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In this paper, we have investigated operator product expansion for thermal correlation function 
of the two scalar currents. Due to breakdown of Lorentz invariance at finite temperature, more 
operators of the same dimension appear in the operator product expansion than those at zero 
temperature. We calculated Wilson coefficients in the short distance expansion and obtain opera- 
tor product expansion for thermal correlation function in terms of quark condensate {ij^ip), gluon 
condensate (Gp^G"^"), quark energy density (uQ^ u) and gluon energy density {uQ^u). 

o 
o 

^ ; I. INTRODUCTION 

! The Shifman, Vainshtein and Zakharov (SVZ) sum rules, proposed about three decades ago [l|, is one of 
the powerful method for investigating the properties of hadrons. This method has been extensively used as an 
C ■ ] efficient tool to study properties of resonances, decay form factors and so on 0. In this approach, hadrons are 
. represented by their interpolating quark currents taken at large virtualities and the correlation function of these 
' currents is investigated in the framework of operator product expansion (OPE). 

. The SVZ sum rules method is extended to the finite temperature in the paper In extending these sum 
O i' rules to finite temperature, two sources of complications arise. One is the interaction of the current with the 
JL. , particles of the medium. The other complication is the breakdown of Lorentz invariance by the choice of the 
reference frame Taking into account both complications, OPE of vector currents at finite temperature first 

" investigated in la] and application of these results to the temperature dependence of the p-meson parameters 
are done in [7|, [p] . Also, modifications of meson parameters due to nuclear medium are widely discussed in the 
] literature [9l.Tl0l|. 

^ , In this paper we studied OPE for thermal correlation function, necessary for the phenomenological investi- 
' gation of scalar meson parameters. We calculated Wilson coefficients in the OPE of scalar currents at finite 
, temperature. 

0^ 
O . 

II. THERMAL CORRELATION FUNCTION IN SHORT DISTANCE EXPANSION 

d 

■ The starting point for SVZ sum rules is the OPE [1] and it gives a general form of the considered quantities 
' in terms of operators. We begin by considering the thermal correlation function. 



>< : T(q) = i / d'xe>"(T(J(x)J{m. (1) 



where J{x) —: q{x)q{x) : is the scalar current and T indicates the time ordered product. The thermal average 
of an operator is defined as follows [11 1 



where H is the QCD Hamiltonian and [3 — 1/T stands for the inverse of the temperature. Traces are carried 
out over any complete set of states. Using Wick's theorem and making some transformations, we get 

T{J{x)J{{))) = trS{x,Q)S{Q,x) - [mx^ : : +2{;i)abX^' : V^.zL'^^f, :] , (3) 

where a and b are spinor indices, and is covariant derivative. The fundamental assumption of Wilson 
expansion is that the product of operators at different points can be expanded as the sum of local operators 
with momentum dependent coefficients in the form: 
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where C„(g^) are called Wilson coefficients and 0„ are a set of local operators. In this expansion, the operators 
are ordered according to their dimension d. The lowest dimension operator with d = is the unit operator 
associated with the perturbative contribution. In the vacuum sum rules low dimension operators composed 
of quark and gluon fields are quark condensate (ipip) and gluon condensate (G^^G'"'^''). At finite temperature 
Lorentz invariance is broken by the choice of a preferred frame of reference and new operators appear in the 
Wilson expansion. To restore Lorentz invariance in thermal field theory, four-vector velocity of the medium u'^ 
is introduced. Using four-vector velocity and quark/gluon fields, we can construct a new set of low dimension 
operators {uQ^u) and (uQ^u) with dimension d = 4. So, we can write thermal correlation function in terms of 
operators up to dimension four: 



T{q) = Cil + C2{-^ij) + G3(G^,G°'^'^) + Ci{uQfu) -f C^iuQ^u) 
In our calculations, we used the following expansion of quark fields '4){x) 



(5) 



(6) 



which takes place in Fock-Schwinger gauge x^A1^{x) — 0. In order to calculate Wilson coefficients in 0], we use 
massless quark field propagator in the background gauge field Q 



S(x,Q) 



1 

4^ 



2 ^ 



B 

48^ 



(7) 



where G^^, = ^A^G"^ is the gluon strength tensor, coefficients a and b have the following forms 



a = 1 



;{u-xf -2lii{-Ax^H^), 



(8) 



6 = 8(M-a;)ln(-4a;V^)- (9) 

In our calculations we get expressions of kind Tap\a — {trGapG\a) , where trace is carry out over color matrices. 
TafjXa is four rank tensor and must be expressed in terms of the metric tensor gap and four- vector velocities Wq . 
The Lorentz covariance at finite temperature allows us to write general structure of this tensor in the following 
form: 

TapXa = Bigapgxa + B2ga\9l3a + Bzgaa9fJ\ + B^gafjUxU^ + B^gaXUpUa + B^gaaUfjUx 

+ BrgpxUaUa + BsgpaUaUx + BggXaUaUp + BioUaUpUxU^ (10) 

where Bi are unknown scalar coefficients. Note that Tapxa- is antisymmetric under the interchange of indices a 
and P, as well as indices A and a: Tapxa — —TpaXa and Tafixa = —TapaX- Using these properties we obtain 
that Bi = i?4 = iJg = i?io = 0, Bq = B-j ^ -B5 = -Bs and B2 = -B3. Therefore Tapxa =< trGapGxa > is 
expressed in terms of two scalar functions of B2 = A and Bq = B. Contracting indices on both sides and using 
gluonic part of energy momentum tensor these coefficients can be expressed with gluon condensate and 
gluon energy densities as follows 

^AGl^G'^n + (11) 



B = Uu^^QlX)- 



(12) 



Using these expressions for coefficients and making some transformations, we obtain the gluonic part of corre- 
lation function 



{trS{x,Q)S{0,x)) 



167r4 



16 aB g 



[u ■ x)hB 2^2 



+ i^Au-xf{uQ^->u) 



(13) 
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Let us consider second term in eq.Q for zero chemical potential case. Expanding the quark bilinear matrix 
in terms of complete set of Dirac matrices and taking into account parity conservation, the most general finite 
temperature decomposition of two operators in spinor space becomes @ 



(ipatpb) = Di6ab + D2{i)ba, 



(14) 



(15) 



Here, Di and E^ are scalar coefficients. By contracting indices on both sides and using fcrmionic part of energy 
momentum tensor O^^,, we get 



(16) 



E2 = liu'^QiX), 



(17) 



Substituting these equalities in eq. ([3]) , the contributions of operators up to dimension four to correlation function 
in the coordinate space can be written as 



{T{j{x)j{m 



1 



27r2a 
1 



-mx'^{ip'ip) — -x^{uQ^u) + -{u • x)'^{u0^u) 

Zi o o 



167r4 

^2 



16 aB {u ■ x)hB g' 



3x4 



8a;' 



2 {G'',.G"n 



29' 



(18) 



Applying the Fourier transformation to eq. (|18p correlation function may be written in the momentum space. 
As can be seen in Fourier transformation, we meet different kind of integrals such as 



d'^x 



(19) 



Up to now, we evaluated expressions in the Minkowski space. After that, in order to calculate integrals we 
proceed to the Euclidean space Using the well known identity 



In a: = — lim 



d 1 



A^o dX (x^)^'' 

the integral in eq. (|19p can be transformed to the following form 



(20) 



I — ^i lim — 

A^o dX 



dx- 



(x2) 



2^l+A ' 



(21) 



where Q is Euclidean momentum and = —q^. Using the standard Fourier transformation and expansion 
formula r(l + e) = 1 — 76 + 0(£^), we get 



87 



1, 9'^ 



(22) 



where 7 is Euler constant. After calculating the other integrals thermal correlation function in the momentum 
space in short distance expansion can be written as 



T{Q) 



QV Q 

2 2 



3 Af4{u-Qf 



327r2Q2\ M 



727r2 
4(u-Q)2 



1 



(5 + 87)7W-(l + 87) 



{n ■ Q? 



1 

4{u-Qf 2 

^Q2 



(23) 
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which reproduces well-known zero temperature result in the T —^ limit |13| . At the end, due to breakdown 
of Lorentz invariance at finite temperature, two new operators appear in the Wilson expansion, in addition to 
the two old (Lorentz invariant) ones, already existing in the vacuum sum rules. Therefore OPE for thermal 
correlator is expressed with four temperature dependent quantities, quark condensate (ipijj): gluon condensate 
(GJJj^G"^'^), quark energy density {uQ^u) and gluon energy density (uQ^u). The thermal average of quark 
condensate (ipip) is known from chiral perturbation theory [13, [iBl 



(VV) = (OIV'V'IO) 



1 T2 



12F2 



(24) 



where n/ is number of quark flavors and F = Q.088GeV. 

The relationship between the trace of energy momentum tensor and the gluon condensate has been studied 
at finite temperature in paper [l6l | 



(e{;) + ^m/((V'V)-(o|^V|o)) 



(25) 



For two massless quarks in the low temperature chira.1 perturbation expansion the trace of the total energy 
momentum tensor Q — UfQ^ + 8^ has following form |i3| 



270 F T 



(26) 



where the pion decay constant has the value of F^r = 0.093GeT^ and the logarithmic scale factor is Ap = 
0.275GeV. The obtained results allow us to investigate the temperature dependence of the a and ao mesons 
parameters. The investigations of hadronic parameters at finite temperature are very important for the inter- 
pretation of heavy ion collision results and for understanding the phenomenological and theoretical aspects of 
thermal QCD. 
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